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Fragmentation functions and their uncertainties are determined for pion, kaon, and proton by 
a global analysis of charged-hadron production data in electron-positron annihilation and by 
the Hessian method for error estimation. It is especially important that the uncertainties of the 
fragmentation functions are estimated in this analysis. The results indicate that the fragmentation 
functions, especially gluon and light-quark fragmentation functions, have large uncertainties at 
small Q^. There are large differences between widely-used functions by KKP (Kniehl, Kramer, and 
Potter) and Kretzer; however, they are compatible with each other and also with our functions if 
the uncertainties are taken into account. We find that determination of the fragmentation functions 
is improved in next-to-leading-order (NLO) analyses for the pion and kaon in comparison with 
leading-order ones. Such a NLO improvement is not obvious in the proton. Since the uncertainties 
are large at small , the uncertainty estimation is very important for analyzing hadron- production 
data at small or pr (Q^, << ) in lepton scattering and hadron-hadron collisions. A code 
is available for general users for calculating obtained fragmentation functions. 



PACS numbers: 13.87.Fh, 13. 66. Be, 13.85.Ni 

I. INTRODUCTION 

In finding any exotic physics signatures beyond the 
standard model and any new hadronic systems in high- 
energy hadron reactions, for example at RHIC (Relativis- 
tic Heavy Ion Collider) and LHC (Large Hadron Col- 
lider), it is necessary to have accurate QCD (Quantum 
Chromodynamics) predictions for cross sections. The 
perturbative QCD is now understood in the next-to- 
leading-order (NLO) level for many reactions and in 
the next-to-next-to-leading order (NNLO) for some pro- 
cesses. Therefore, the crucial part is to understand non- 
perturbative aspects, namely parton distribution func- 
tions (PDFs) and fragmentation functions. 

Fragmentation functions arc used in high-energy re- 
action processes with hadron production. They include 
hadron-production processes in electron-positron annihi- 
lation, electron, muon, and neutrino scattering from pro- 
ton or nucleus, proton-proton collisions, and heavy-ion 
collisions [l|, Q. Such reactions are becoming increas- 
ingly important in hadron physics for investigating the 
origin of the proton spin and quark-hadron matters in 
heavy-ion reactions. 

The fragmentation functions are related to a nonper- 
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turbative aspect of QCD, so that they cannot be precisely 
calculated by theoretical methods at this stage. The sit- 
uation is similar to the determination of the PDFs, where 
high-energy experimental data are used for their deter- 
mination instead of theoretical calculations. There are 
already several studies on such analyses 0] . Widely used 
parametrizations were obtained by Kniehl, Kramer, and 
Potter (KKP) [1] and also by Kretzer [Hj. An updated 
version was reported recently by Albino, Kniehl, and 
Kramer (AKK) The fragmentation functions were 
determined by analyzing the hadron-production data in 
the e^e~ annihilation. The analysis was done in leading 
order (LO) and NLO of nmning coupling constant a^. 

Despite the importance of the fragmentation functions, 
it is unfortunate that uncertainties of the fragmentation 
functions have not been estimated yet. The uncertainties 
have been investigated extensively for unpolarized PDFs 
6], polarized PDFs and nuclear PDFs In order 
to extract any information from the hadron-production 
processes, reliability regions of the fragmentation func- 
tions should be shown in comparison with actual data. 
In particular, it is known that there are large differences 
between the functions of KKP and Kretzer, for example 
between their gluon functions. An error analysis should 
be done for the fragmentation functions. 

In this paper, we analyze the e^e~ data for obtain- 
ing the fragmentation functions and their uncertainties 
in both LO and NLO. In particular, this analysis is the 
first attempt for showing the uncertainties of the frag- 
mentation functions . Therefore, our works focus on 

• determination of the fragmentation functions and 



their uncertainties for pion, kaon, and proton in 
both LO and NLO, 

• roles of NLO terms on the fragmentation-function 
determination, namely on their uncertainties, 

• comparison with other analysis results by consider- 
ing the uncertainties. 

The functions are defined at an initial point in terms 
of a number of parameters, which are determined so as 
to explain the e"*" -I- e~ h + X data, where h = tt^, 
K^, and p/p. The uncertainties are calculated by the 
Hessian method, which has been used for obtaining the 
PDF uncertainties in Refs. [i,0,[il- 

This paper is organized as follows. A formalism is dis- 
cussed in Sec. |ll]for hadron-production cross sections in 
the e'^e^ annihilation with the fragmentation functions. 
In Sec. mil our analysis method is described for deter- 
mining the fragmentation functions. Analysis results are 
explained in Sec. IIVI The results are summarized in Sec. 
Ivl In Appendix |Al we explain how to use the obtained 
functions to other hadrons, 7r°, K'^ , , n, and n. 



II. FORMALISM 



where the perturbative correction is included up to the 
NLO. In the LO analysis, the second term is removed. 
The electroweak cross section for producing a qq pair is 
given by [ll| 
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where the terms with pi(s) and P2(s) come from j-Z 
interference and Z processes, and they are given by 
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Here, a is the fine structure constant in Quantum Elec- 
trodynamics (QED), eq is a quark charge, and 6w is the 
weak-mixing angle. The mass and width of Z are denoted 
by Mz and F^, respectively. The vector and axial- vector 
couplings Cy and c;^ of a fermion / are expressed by the 
third component of the weak isospin and the fermion 
charge e/: 



A. Hadron production in e+e annihilations 
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The cross section for the reaction e+ + ^ ft. + AT is 
described by two step processes Q. The first part is to 
create a quark-antiquark pair by the reaction e~^e~ — *■ qq, 
and higher-order corrections such as e~^e~ qqg are also 
taken into account in the NLO analysis. The second part 
is to create a hadron h from quark (g), antiquark [q), or 
gluon ((7), and this process is called fragmentation. 

The fragmentation function is defined by the hadron- 
production cross section and the total hadronic cross sec- 
tion [ii: 
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where is the virtual photon or Z momentum squared 
in e~'"e~ —> Z and it is given hy — s with the 
center-of-mass energy y^. The variable z is defined by 
the energy fraction: 
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where Eh and \/s/2 are the hadron and beam ener- 
gies, respectively. Namely, it indicates the hadron energy 
scaled to the beam energy. 

The total cross section is described by the qq-paii cre- 



ation processes, e"'"e 
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with higher-order corrections: 
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Actual expressions in terms of sin 9w are 
-i for the electron, 

for up, charm, and top quarks, 
^ + 1 sin^ 9w and 



and Cy 



and bottom quarks Q. 



i for down, strange, 



B. Fragmentation functions and their Q evolution 

The fragmentation process occurs from primary 
quarks, antiquarks, and gluons, so that F^{z,Q'^) is ex- 
pressed by the sum of their contributions ^lOj : 



(7) 



where D^{z,Q^) is a fragmentation function of the 
hadron h from a parton i (= u, d, s, ■ ■ •, g), C'i{z,as) 
is a coefficient function, and the convolution integral 
is defined by 



y \y 
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The function D^{z, Q^) indicates the probability to find 
the hadron h from a parton i with the energy fraction z. 
The coefficient functions are calculated in perturbative 
QCD, and the NLO results are listed, for example, in 
Refs. [3, S [Hi for the modified minimal subtraction 
(MS) scheme. The cross section is split into longitudinal 
and transverse components for the virtual 7 or Z, so 
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that the coefficient functions are also expressed by these 
components. 

evolution for the fragmentation functions is cal- 
culated in the same way as the one for the PDFs by 
using the timelike DGLAP (Dokshitzer-Gribov-Lipatov- 
AltarcUi-Parisi) evolution equations. For example, the 
flavor-singlet evolution is given by 



(9) 



where Dg{z,Q'^) denotes the singlet function 
D^si^.Q^) = E,[D'i{z,Q') + Dl{z,Q^)], and Nf 
is the number of quark flavors. One should be careful 
that the off-diagonal elements Pgq{z) and Pqg{z) are 
interchanged in the splitting-function matrix from the 
PDF case. The splitting functions are the same in the 
LO evolution of the PDFs; however, they are different 
in the NLO. Explicit forms of the splitting functions are 
provided in Refs. [lol.[l^. 

The evolution equations are essentially the same as 
the PDF case, so that the same numerical method can 
be applied for obtaining a solution. The equations are 
solved by direct integrations in the z space as explained 
in Ref. [13] . A slightly modified numerical approach with 
the Simpson's integration method is used in this analysis, 
and its evolution results are independently checked by 
using a Gauss-Legendre quadrature. Actual evolutions 
are calculated for the functions — D'^ ± Do, 
and Z3g. Then, the quark and antiquark functions are 
obtained by their combinations: D'^ — {D\ + D'^_)/2 
andi?,^ = m\-i^„^)/2. 



III. ANALYSIS METHOD 

A. Parametrization 

The fragmentation functions are expressed in terms of 
a number of parameters at the initial (= Q^) in the 
same way as the PDF analysis [1, H [1| ■ Since they should 
vanish at z=l, a simple polynomial form is taken: 



D^{z,Ql)^N^z-Hl~zfl 



(10) 



where , a'', and (3^ are parameters to be determined 
by a analysis of e+e~ hX data. This kind of poly- 
nomial form has been assumed in the parametrization 
studies of the fragmentation functions [j, Q . The scale 
Qq is taken at mass thresholds and for charm 
and bottom functions, where rric and mi, are charm- and 
bottom-quark masses, respectively. For the hadron h, we 
have pious (7r+-|~7r~), kaons {K~^+K~), and proton/anti- 
proton (j> + p). 



From actual analysis trials, we found that it is more 
appropriate to take the second moment A//' as a parame- 
ter instead of N^'-. This is because there exists an energy 
sum rule 
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for the function D^{z, Q^), so that M/' should not exceed 



< 1. 



(12) 



This constraint should be imposed in a global analysis. 
Another advantage is that the physical meaning of 
is clear. It is the energy fraction for the hadron h which 
is created from the parton i. These parameters iV/* and 
are related with each other by 



TV'' = 



(13) 



where B(a^ -I- 2, f]'^ + 1) is the beta function. 

A general principle of our parametrization is to use 
a common function for favored fragmentation functions 
from up and down quarks. A separate function is used 
for a favored one from a strange quark by considering the 
mass difference. Different functions are assigned for dis- 
favored ones. Here, the favored means the fragmentation 
from a quark which exits in the hadron h as a constituent 
in the naive SU(6) quark model. The disfavored means 
the fragmentation from a sea quark. A flavor symmet- 
ric form is assumed for the fragmentation functions from 
light sea-quarks (up, down, and strange sea-quarks). Al- 
though it is known that light sea-quark distributions are 
not flavor symmetric in the unpolarized PDFs [l3| , there 
is no data to distinguish among the fragmentation func- 
tions from different light sea-quarks. 

The actual parametrization forms are shown for the 
pion, kaon, and proton in the following. 

(1) Pion (7r+) 

Considering the constituent quark composition 
TT^{ud), we take the same favored fragmentation 
functions for tt"*" from u and d quarks: 

D:-" (z, Ql) = Df [z, Ql) = Nf z<^ (1 - zr^^ . (14) 

The pion productions from u, d, s, and s are disfavored 
processes, and they are considered the same at the initial 
scale: 

Df{z,Ql)=Df{z,Ql) 



Di {z.qD^di (z,g2) 

Nfz<\l-zfl\ 



(15) 



In addition, a fragmentation function from gluon is given 
by 



Df{z,Ql)^Nfz< (l-z)'^r 



(16) 
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These functions are provided at the initial scale Qq with 
the parameters. 

Different functions are assigned for productions from 
heavy quarks because of mass differences: 

Df{z,ml) = Df{z,ml) = Nf z^'^\l - z f^\ (17) 
Dfiz,ml) = Df{z,ml) = N^\<\l - zf^\ (18) 

Thresholds for heavy quarks are — rri^ and in cal- 
culating evolutions and the running coupling constant 
as{Q^)- However, the thresholds for the cross section are 
taken = Ami and Ami f^i 5]. 

At different (7^ Qg); one should note that some 
equal relations do not hold in the NLO because of a g ^ g 
or q q splitting [l3| : 

Df (z, Q2) = Df (z, Q2) ^ Df (z, Q'), (19) 

although i?;+(z,Q2) = Dfiz,Q^), Df{z,Q^) = 
Df{z,Q^), Dfiz,Q^) = Df{z,Q^), and 
D^^ {z, Q2) ^ of (z, Q2) are still valid. 

(2) Kaon {K+) 

Parameters are assigned to A'+ fragmentation func- 
tions in the same way by considering the constituent 
quark composition K~^{us). The only difference is that 
the anti-strange function is taken in a different form from 
the up-quark one: 

Dr{z,Ql)=Nrz-^Ul-zf^\ (20) 

Driz,Qi)-Nrz-^\i-zr^\ (21) 

(z, Ql) = of {z,Ql)= Df {z,Ql) 

^Dr{z,Ql)-Nrz-'^\l-zf^\ (22) 

or iz, ml) ^ dT (z, ml) = iVf (1 - z)^!^" , 

(23) 

Df [z, ml) = Df (z, ml) = Nr z<" (1 - z)^^" , 

(24) 

Df{z,Ql)=Nfz<\l-z)^r. (25) 

At general (^ Qq), there are differences in the NLO 
as explained in the pion case: 

Driz,Q') ^ Dfiz^Q^) ^ Dfiz,Q^). (26) 

The relations Df{z,Q^) = Df{z,Q'^), Df{z,Q^) = 
Df{z,Q^), and Df{z,Q^) = Df{z,Q^) are stiU 
valid in the evolution. 



(3) Proton (p) 

Proton fragmentation functions are also parametrized 
in the same way by considering the constituent quark 
composition p (uud): 

Dl{z,Ql) = 2D2z,Ql) = NPz<{l - z)''£, (27) 
DUz,Ql) = DP^{z,Ql) ^ DPiz,Ql) 

= DUz,Ql) = Niz<{l-z)f''-, (28) 

DP{z,ml) = DUz,ml) = iV^'z<(l - z)'^?, (29) 

Dl{z,ml) = Dl{z,ml) = Nlz<{l - z)''s (30) 

Dl{z,Ql)^Nlz<{l^z)Pl. (31) 

The major difference from the mesons is the factor of two 
in £'S(z, Qq) — 2D^{z,Qq) which is suggested simply by 
considering valence-quark structure with a flavor sym- 
metry In order to produce a baryon from a quark, 
two qq pairs need to be created. If the initial quark is 
up quark, two creation process, {uu){dd) and (dd){uu), 
should contribute to the proton formation, whereas only 
the (uu) (uu) process contributes if the down quark is in 
the initial state. It leads to the factor of two although it 
is a naive counting estimate. This factor is also the same 
in a spectator di-quark model in the flavor symmetric 
case [l5| . 

At different Q^, one should be careful about the fol- 
lowing inequalities in the NLO: 

DS(z,Q2)^2D^(z,g^), (32) 
DUz, g2) ^ D|(z, Q2) ^ DP(^z, Q\ (33) 

The equalities for strange, charm, and bottom quarks, 
i^f(z,Q2) = 7?f(z,Q2), i?P(z,Q2) = i5f(z,Q2), and 

Dl{z,Q'^) = D?{z,Q'^), are not changed by the evo- 
lution. 

(4) TT^, , and p 

In the constituent quark model, 7r~, K~ , and p are ex- 
pressed Tr~ (ud), K~{us), and p{uud). Using the charge 
symmetry, we relate the fragmentation functions of 7r~ 
to the ones of 7r+: 

Dl^^-^nz,Q^) = Df^''"^^{z,Q^). (34) 

The gluonic functions are the same: 

Dl^^-^f>{z,Q^) = Dl^^^^^^{z,Q^). (35) 

The functions for 7r°, , , n, and n are also obtained 
by using the charge symmetry, and it is explained in Ap- 
pendix 1^ 

After all, the following parameters are used for the pion 
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kaon, and proton in our global analysis: 



pion: af, Pl\ Alf, af, (3f, Alf, 

af, pf, Mf,al\ pf, Mf, af, pf, Alf, 
kaon- a^^ 0^^ M^^ a^^ 3^^ M^^ a^^ 3^^ M^^ 

proton: P^, M^, al, Pi, Ml, 



B. Experimental data 

The fragmentation functions are determined by the 
charged- hadron production data of e+ + e~ ^ /i^ + X. 
Used data are those from the measurements of TASSO 
[i1[I1[i|L TPC dljiRS H, TOPAZ J21LSLD [H, 
ALEPH (H, OPAL [il, and DELPHI [2^j2^. There 
are data from MARK-II ^ and JADE [2l| collabora- 
tions; however, they are not used in our analysis because 
of the kinematical condition z > 0.1, which is explained 
in the following. References, center-of-mass (cm.) ener- 
gies, and numbers of the data are listed for used data sets 
in Tables mini and |TTI] ;29] . The ALEPH, DELPHI, and 
OPAL data are taken at CERN (European Organization 
for Nuclear Research), TASSO at DESY (German Elec- 
tron Synchrotron), TOPAZ at KEK (High Energy Accel- 
erator Research Organization), and HRS, SLD, and TPC 
at SLAG (Stanford Linear Accelerator Genter). In the 
DELPHI and SLD measurements, the light-quark and 
heavy-quark events are separated. 

Since the perturbative QGD is applied in the evo- 
lution calculations, the data with > 1 GeV^ are used 
in our analysis. Furthermore, small-z data are excluded 
because of soft-gluon emission. The minimum z value is 
Zmin =0.1 for the data at ^/3 < Mz and Zmin — 0.05 
for the data at ^/s = Mz- Resummation effects of soft- 
gluon logarithms need to be clarified in order to include 
the small-z data into the analysis [l^l . We do not include 
three-jet and unidentified-hadron data [1, [3l|. There are 
ambiguities for extracting the gluon fragmentation func- 
tion from three-jet cross sections 5]. In oder to describe 
the unidentified-hadron cross sections, all the relevant 
hadrons need to be included in the analysis, whereas only 
the pions, kaons, proton, and anti-proton are taken into 
account in our analysis. 

The kinematical range is shown in Fig. [l]by the vari- 
ables z and Q{= \fQ^ — \fs) for the pion data. Here, 
the data are also shown in the small z region (z < 0.1 
or z < 0.05), where they are not included in the actual 
analysis. The range of Q is wide from the low-energy 
TASSO data to high-energy SLD, ALEPH, OPAL, and 
DELPHI data at Q = Mz- Many data are collected in 
the small- 2: region [z < 0.4), whereas the large-z data are 
scarce because cross sections are very small. 



TABLE I: Experiments, references, center-of-mass energies, 
and numbers of data points are listed for used data sets of 



+ e' 



+ X 



experiment 


rcf. 




# of data 


TASSO 


[16, 17, 18112,14,22,30,34,44 


29 


TPC 


[19] 


29 


18 


HRS 


[20] 


29 
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58 
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SLD 


[22] 


91.28 


29 


SLD (u,d,s quark) 


[22] 


91.28 


29 


SLD (c quark) 


[22] 


91.28 


29 


SLD (b quark) 


[22] 


91.28 


29 


ALEPH 


[23] 


91.2 


22 


OPAL 
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91.2 


22 


DELPHI 


[25] 


91.2 


17 


DELPHI (u,d,s quark) 


[25] 


91.2 


17 


DELPHI (b quark) 


\M 


91.2 


17 


total 






264 


TABLE II: Experimental information 


is listed for the used 


data of e+ + e K 








experiment 


ref. 




# of data 



TASSO 

TPC 

HRS 

TOPAZ 

SLD 

SLD (u,d,s quark) 
SLD (c quark) 
SLD (b quark) 
ALEPH 
OPAL 
DELPHI 

DELPHI (u,d,s quark) 
DELPHI (b quark) 



[16, rr,18] 12,14,22,30,34 



[19] 
20] 

22] 
22] 
22] 
32] 
23] 

m 

25, 26] 
[25] 
[25] 



29 

29 

58 
91.28 
91.28 
91.28 
91.28 
91.2 
91.2 
91.2 
91.2 
91.2 



18 
17 
3 
3 
29 
29 
29 
28 
18 
10 
27 
17 
17 



total 



245 



TABLE III: Experimental information is listed for the used 
data of e+ +e- -^p/p + X [2^. 
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24 
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[19] 
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17 
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[22] 
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SLD (u,d,s quark) 


[22] 
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SLD (c quark) 
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[22] 
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[23] 
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18 
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10 
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[25, 26] 


91.2 


23 


DELPHI (u,d,s quark) 


[25] 


91.2 


17 


DELPHI (b quark) 


[25] 


91.2 


17 


total 






249 
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FIG. 1: (Color online) Kinematical range is shown by z and 
Q(= values for the pion data. 



C. X analysis 

The total is calculated in comparison with the data 
for the fragmentation function F'^{z,Q'^) of the hadron 
h in Eq. ([T]). The scale is equal to the cm. en- 
ergy squared s in Tables HI [TTl and IIIII The initial scale, 
where the fragmentation functions for up, down, and 
strange quarks are defined, is taken Q^—l GeV^. Since 
the e+e" data are measured at large as shown in 
the tables, we could have taken a larger Qo i^i the range 
where perturbative QCD can be applied undoubtedly 
without worrying about higher-twist and higher-order 
corrections. However, the fragmentation functions are 
practically used also in a relatively small region, for 
example, in hadron-production process in deep inelastic 
lepton scattering [l| and proton-proton collisions • For 
applications of obtained functions to such processes, it is 
better to supply a code for calculating the functions even 
at small (~ 1 GeV^). Considering this situation, we 
decided to take Qo=l GeV^. The heavy-quark functions 
are defined at the mass thresholds. The quark masses are 
?7ic=1.43 GeV and irib—A.S GeV as given in Ref. [s^ . 

The theoretical functions should be obtained at the 
same experimental points for calculating x^- The 
evolution is calculated by the timelike DGLAP equations 
as explained in Sec. IIIBI Then, the evolved functions 
are convoluted with the coefficient functions in Eq. ([7]). 
These calculations depend on the choice of the scale pa- 
rameter A. It is taken A^'^^=0.220 GeV or a|^^^o=0.323 
GeV for four flavors [s^l. The MS scheme is used in the 
NLO. Below the charm threshold and above the bottom 
one, these values are converted to corresponding ones for 
three and five flavors so that the coupling constant is 
continuous at the thresholds. 

The total is then calculated by 



E 



data 
3 



(36) 



imental point. The experimental errors are calcu- 
lated from systematic and statistical errors by (o-j'''*")^ = 
(cr|^*)^ _l_ ^^staf^2_ rpj^g assigned parameters are deter- 
mined so as to obtain the minimum . The optimization 
of the functions is done by the CERN subroutine MINUIT 



D. Uncertainties of fragmentation functions 

Uncertainties of the PDFs were estimated in Refs. 
[1, 0, Q • Two methods have been developed: the Hes- 
sian and Lagrange-multiplier methods. Although tech- 
nical details are described in these references, outline of 
the Hessian method is explained because it is used in our 
analysis and uncertainties have not been investigated in 
the fragmentations functions. 

The parameters are denoted 2, • • •, N), where 

N is the number of the parameters. We expand 
around the minimum point ^: 

Ax'(0 = X'ii + SO - x\i) = E ii^M^Si, , (37) 

where only the leading quadratic term is kept, and the 
Hessian Hij is the second derivative matrix. The confi- 
dence region is given in the parameter space by supplying 
a value of Ax^- It is known that the confidence level is 
68% for Ax^=l if the number of the parameters is one 
[N — \). The Ax^ value needs to be changed in a gen- 
eral case of ^ 1. Assuming correspondence between 
the confidence level of a normal distribution in multi- 
parameter space and the one of a x^ distribution with N 
degree of freedom, we have the confidence level P: 



P 



1 



2 r(A^/2) 



exp 



dS , (38) 



where r(A^/2) is the Gamma function. The value of Ax'^ 
is taken so that the confidence level becomes the one-cr- 
error range, namely P = 0.6826 0, Q. The details of 
this Ax^ choice are explained in Ref. jsj. The Ax^ 
value is numerically calculated by using Eq. ([55]) . For 
example, Ax^ = 15.94 is obtained for A^ — 14. The 
Hessian matrix is obtained by running the subroutine 
MINUIT. From this Hessian, Ax^, and derivatives of the 
fragmentation functions with respect to the parameters, 
the uncertainties are calculated by 



[dD':iz)r = Ax^ 



3,k 



do 



H 



Jk 



\ do 



(39) 



where Fj and are experimental and theoreti- 

cal values of F''-{z,Q'^), respectively, at the same exper- 
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IV. 



RESULTS 



TABLE IV: Parameters determined for the pion. 



Analysis results are explained. In Sec. IIV A| opti- 
mized parameters are shown, and contributions from 
data sets are listed. Then, fit results are compared with 
experimental data. In Sec. IIVB| the obtained fragmen- 
tation functions and their uncertainties are shown. They 
are compared with other parametrization results in Sec. 
HVCl 



A. Comparison with experimental data 

Obtained parameters in the LO and NLO are listed 
in Tables IIVI IVl and IVII for the pion, kaon, and proton, 
respectively. In these analyses, it was very difficult to 
determine the gluon functions, so that we decided to fix 
some parameters. In trial analyses, we found that the 
gluon function for the pion {zD^^) tends to be peaked 
at z ~ 0.2. The parameter (3g controls its functional 
behavior at large z. Since the gluon function becomes 
very small at large z, the value of Pg does not affect 
the value to a significant amount if the function is 
peaked at small z. Therefore, the parameter Pg is fixed 
as (3g=S in the pion analyses. If this value is taken much 
more than eight, the energy sum rule in Eq. (jlip is badly 
violated in the gluon part. If it much less than eight, it is 
violated in the up-quark part. The value Pg=S> is chosen 
to compromise these two issues. In the analysis of Kretzer 
[1| , one of the gluon parameters is also fixed because the 
second moment is equal to the average of the moments 
of up- and down-quark fragmentation functions. 

In the kaon analyses, the gluon function tends to be 
peaked at large z. Therefore, ag becomes the parameter 
which does not affect the total x^ instead of (3g. It is 
chosen = 5 so that the moments for gluon and up 
quark do not become too large to affect the sum rule to 
a significant extent. 

It was also very difficult to determine the gluon dis- 
tribution for the proton. It indicates that the data 
are not sensitive to the gluon function at this stage 
even in the NLO analysis. Therefore, the moment of 
the gluon function is fixed at the average of the mo- 
ments for the favored and disfavored functions Q : Mg = 
[{Mu + Md)/2 + Ms]/2. Furthermore, the function tends 
to be peaked at large z^ so that the parameter Ug is fixed 
at ag = 5 in the same way as the kaon case. The average 
of the moments is motived by the following consideration. 
The process for producing a proton from an initial up or 
down quark should contain two quark-pair productions 
(uu and dd pairs or two uu pairs), so that it is propor- 
tional to with the strong coupling constant g. In the 
same way, the process from an initial u quark is propor- 
tional to , and the one from an initial gluon is to g^ . 
Therefore, the gluon moment is, roughly speaking, given 
by the average of favored {u and d) and disfavored (u) 
moments. Such averages are also taken in Ref. @. 



function 



M 



(LO) 

of 
of 
of 



0.546 ± 0.085 
0.250 ± 0.068 
0.305 ± 0.046 
0.302 ± 0.023 
0.115 ± 0.111 



-1.100 ± 0.183 
-0.500 ± 0.301 
-1.007 ± 0.123 
-1.176 ± 0.045 
1.405 ± 0.897 



1.282 ± 0.140 
5.197 ± 0.576 
3.918 ± 0.236 
5.805 ± 0.188 
8.0 (fixed) 



(NLO) 
Df 

of 



0.401 ± 0.052 
0.094 ± 0.029 
0.178 ± 0.018 
0.236 ± 0.009 
0.238 ± 0.029 



-0.963 ± 0.177 
0.718 ± 0.466 
-0.845 ± 0.108 
-1.219 ± 0.042 
1.943 ± 0.399 



1.370 ± 0.144 
6.266 ± 0.808 
3.868 ± 0.323 
5.668 ± 0.219 
8.0 (fixed) 



TABLE V: Parameters determined for the kaon 



function 


M 


a 


P 


(LO) 


0.0922 ± 0.0419 


0.588 ± 1.605 


1.632 ± 0.812 


Dt 


0.1651 ± 0.0962 


2.190 ± 2.871 


2.829 ± 1.143 


Dt 


0.0638 ± 0.0363 


0.565 ± 0.460 


7.093 ± 3.383 


or 


0.0919 ± 0.0055 


0.230 ± 0.157 


4.549 ± 0.388 


or 


0.0696 ± 0.0027 


0.017 ± 0.110 


8.808 ± 0.534 


of 


0.0319 ± 0.0147 


5.0 (fixed) 


0.247 ± 0.162 


(NLO) 

jjK+ 


0.0740 ± 0.0268 


-0.630 ± 0.629 


1.310 ± 0.772 


^ s 


0.0878 ± 0.0506 


2.000 ± 2.913 


2.800 ± 1.313 




0.0255 ± 0.0173 


1.012 ± 0.939 


8.000 ± 3.715 


or 


0.0583 ± 0.0052 


0.527 ± 0.228 


5.866 ± 0.636 


or 


0.0522 ± 0.0024 


0.247 ± 0.126 


11.212 ± 0.721 




0.0705 ± 0.0099 


5.0 (fixed) 


0.810 ± 0.239 



TABLE VL Parameters determined for the proton. 



function 



M 



P 



(LO) 
Dl 



0.0839 ± 0.0125 
0.0158 ± 0.0058 
0.0241 ± 0.0015 
0.0180 ± 0.0006 

(Mt.+Mrf)/2+A/fl 



-0.814 ± 0.200 
0.866 ± 0.677 
0.683 ± 0.359 
0.071 ± 0.178 
5.0 (fixed) 



1.628 ± 0.324 
5.078 ± 1.400 
7.375 ± 1.250 
8.802 ± 0.839 
2.927 ± 1.117 



(NLO) 

Dl 

Dl 

Dl 

Dl 

Dl 



0.0732 ± 0.0113 
0.0084 ± 0.0057 
0.0184 ± 0.0017 
0.0155 ± 0.0007 

(M^+Mrf)/2+Mfl 



-1.022 ± 0.219 1.434 ± 0.268 

1.779 ± 1.422 4.763 ± 1.882 

0.407 ± 0.373 6.784 ± 1.555 

-0.203 ± 0.165 8.209 ± 0.950 

5.0 (fixed) 4.900 ± 2.046 
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TABLE VII: Each contribution in the pion analysis. 



6xp6riiiiGnt 


# of data 








9Q 


^9 1 


OL.iJ 


TPP 




oo.o 




HRS 


2 


1.1 


2.0 


TOPAZ 


4 




9 fi 


SLD (all) 


29 


11.3 


10.6 


SLD (u,d,s) 


29 


46.0 


36.4 


SLD (c) 


29 


24.4 


26.1 


SLD (b) 


29 


71.2 


66.4 


ALEPH 


22 


22.8 


24.0 


OPAL 


22 


45.4 


45.8 


DELPHI (all) 


17 


48.3 


48.6 


DELPHI (u,d,s^ 


) 17 


29.6 


31.1 


DELPHI (b) 


17 


64.9 


60.8 


total 


264 


453.2 


433.5 


(/d.o.f.) 




(1.81) 


(1.73) 


TABLE VIII: 


Each contribution in the kaon analysis. 


pynpriiTiPnf 


# of data 


(LO) 


x' (NLO) 




18 


26.8 


25.0 




17 


15.2 


15.2 


HRS 


3 


0.6 


0.4 


TOPAZ 


3 


0.5 


0.8 


SLD (all) 


29 


14.9 


12.3 


SLD (u,d,s) 


29 


58.8 


57.2 


SLD (c) 


29 


33.6 


32.4 


SLD (b) 


28 


127.7 


88.7 


ALEPH 


18 


10.3 


12.8 


OPAL 


10 


10.6 


11.5 


DELPHI (all) 


27 


14.8 


15.2 


DELPHI (u,d,s 


) 17 


22.5 


22.1 


DELPHI (b) 


17 


11.7 


11.7 


total 


245 


348.2 


305.1 


(/d.o.f.) 




(1.53) 


(1.34) 



TABLE IX: Each contribution in the proton/anti-proton 
analysis. 



experiment 


# of data 


X' (LO) 


x' (NLO) 


TASSO 


24 


34.9 


33.6 


TPC 


17 


22.6 


23.5 


HRS 


4 


19.5 


16.3 


TOPAZ 


3 


3.6 


3.2 


SLD (all) 


30 


13.0 


12.6 


SLD (u,d,s) 


30 


62.9 


56.6 


SLD (c) 


30 


46.3 


47.1 


SLD (b) 


26 


36.3 


37.9 


ALEPH 


18 


15.7 


15.8 


OPAL 


10 


112.1 


110.6 


DELPHI (aU) 


23 


5.6 


6.4 


DELPHI (u,d,s) 


17 


2.1 


2.4 


DELPHI (b) 


17 


16.5 


17.1 


total 


249 


391.2 


383.2 


(/d.o.f.) 




(1.66) 


(1.62) 



The difficulty in determining the gluon function is re- 
flected in large errors of the gluon parameters as shown 
in Tables HvllVl and ED The hght-quark (u, u, s) func- 
tions also have large errors which are as large as the gluon 
ones. The values of the second moments M/* indicate 
that a large fraction for the final hadrons comes from the 
pions. Kaon and proton contributions are rather small. 
The moments also indicate that the favored functions 
are generally larger than the disfavored ones as expected. 
Adding these moments, we find that the total moments 
are generally within the energy sum rule in Eq. pip . 
However, the sums slightly exceed one for the gluon and 
up quark. Considering the errors in the moments, we did 
not strictly impose the sum-rule condition. We found 
a general tendency that the moment errors for the pion 
become smaller in the NLO in comparison with the LO, 
whereas the NLO errors are as large as the LO ones for 
the kaon and proton. Although many used data are taken 
at large as shown in Fig. [1] the pion data are sen- 
sitive to the NLO corrections in the evolution and 
coefficient functions. 

Each contribution is listed in Tables IVIIl IVIIIl and 
IIXI for the pion, kaon, and proton/anti-proton, respec- 
tively. The values indicate that the pion data of HRS, 
TOPAZ, SLD (all; u,d,s; c), and ALEPH are explained 
by our parametrization, whereas TASSO, TPC, SLD (5), 
OPAL, and DELPHI (all; u,d,s; b) data are not so well re- 
produced. In the kaon analysis, the situation is slightly 
different. The data of TPC, HRS, TOPAZ, SLD (all; 
c), ALEPH, OPAL, and DELPHI (all; u,d,s; b) are now 
well explained by our fits, whereas the data of TASSO 
and SLD (u,(i,s; b) deviate. In the proton analysis, the 
data of TPC, TOPAZ, SLD(all), ALEPH, and DELPHI 
(all; u,d,s; b) are well explained; however, the data of 
TASSO, HRS, SLD {u,d,s; c; 6), and OPAL are not re- 
produced. Some data sets deviate significantly from our 
parametrization, namely from other data sets. For ex- 
ample, the x^ value is more than hundred for ten data 
points in the OPAL proton data. 

In Fig. [21 our NLO parametrization results are com- 
pared with the data at Q = Mz without separation on 
initial partons by the ALEPH, DELPHI, OPAL, and SLD 
collaborations. The figure shows the comparison with the 
pion data for F'^ (z, Q'^) in Eq. ([1]). Here, tt^ indicates 
charged pions: tt^ = 7r"'"-|-7r~ . The one-cr range of the un- 
certainty band is shown by the shaded area. The theoret- 
ical curve and its uncertainty are calculated at Q = Mz- 
Most of the data are explained by the parametrization 
and they are within the uncertainty band. However, some 
DELPHI data are outside the band, which leads to the 
large x^ contribution in Table IVTIl 

Next, the NLO results are compared with each data 
set for the pion in Figs. [3] and 21 where the data 
for the initial-quark separation are also shown. The 
figure indicates the rational difference between the 
data and the theoretical parametrization: [F'^ (data) — 
F^ (theory)] /F'^ (theory) at the same point with 
the data. Most of the data are compatible with the NLO 
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FIG. 2: (Color online) Comparison of our NLO results with 
pion-production data at Q = Mz without separation on ini- 
tial partons by the ALEPH, DELPHI, OPAL, and SLD col- 
laborations. 



parametrization; however, some DELPHI data are out- 
side the range of our theoretical estimations. There may 
be inconsistency among the data sets. However, we did 
not remove the DELPHI data from the analysis because 
the obtained functions and the total did not change 
significantly even if they are excluded from the data set. 
The heavy-quark data also deviate slightly from the the- 
oretical calculations, which result in the relatively large 
contributions. 

The kaon fit results are shown in Figs. [5] and El The 
kaon data are reproduced well as shown in these figures. 
There is no data which indicates significant deviation 
from the uncertainty regions, which is reflected in the 
smaller (x^/d.o.f.=1.34 in NLO) than the value in 
the pion analysis (x^/d.o.f.— 1.73). However, the bottom- 
quark data by the SLD are not well explained in the re- 
gion, 0.4 < z < 0.5. 

The proton results are shown in Figs. [7] and [51 Most 
of the data are within the uncertainty bands. However, 
there are serious deviations for the OPAL data, which 
results in the huge value (x^ = 110.6 in NLO) for only 
ten data points. Since only the OPAL data are different 
from the other ones by the ALEPH, DELPHI, and SLD 
collaborations, there is an inconsistency problem for the 
OPAL measurements. 

These figures indicate that our fits arc successful. The 
uncertainty bands become large at small- and large-z re- 
gions, so that experimental measurements and theoret- 
ical studies are needed to determine the fragmentation 
functions in the wide-z range. 



TASSO Q=12GeV _ 14GeV t 22GeV 

T 




+1 

I 

'a 
■o 



■ SLD H OPAL DELPHI 
^> ALEPH 5 

2 ^ _™,«n»»inTn 




^^5- -g:Aag^iijlluj|ili]]^W|| 

Q = Mz 






FIG. 3: (Color online) Comparison with charged-pion pro- 
duction data by the TASSO, TPC, HRS, ALEPH, DELPHI, 
OPAL, and SLD collaborations. The rational differences be- 
tween the data and theoretical calculations are shown as a 
function of z. The average scale Q=16 GeV is taken for the- 
oretical calculations in the top figure with the TASSO data. 




FIG. 4: (Color online) Comparison with charged-pion produc- 
tion data by the TASSO, TOPAZ, DELPHI, and SLD collab- 
orations. The scale is Q—39 GeV for theoretical calculations 
in the top figure. 
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FIG. 5: (Color online) Comparison with charged-kaon pro- 
duction data by the TASSO, TPC, HRS, ALEPH, DELPHI, 
OPAL, and SLD collaborations. The rational differences be- 
tween the data and theoretical calculations are shown. The 
scale is Q—16 GeV for theoretical calculations in the top fig- 
ure. 



FIG. 7: (Color online) Comparison with proton/anti- proton 
production data by the TASSO, TPC, HRS, ALEPH, DEL- 
PHI, OPAL, and SLD collaborations. The rational differences 
between the data and theoretical calculations are shown. The 
scale is Q=16 GeV for theoretical calculations in the top fig- 
ure. 




ft, 
I 




FIG. 6: (Color online) Comparison with charged-kaon pro- 
duction data by the TASSO, TOPAZ, DELPHI, and SLD 
collaborations. The scale is Q=34 GeV for theoretical calcu- 
lations in the top figure. 



FIG. 8: (Color online) Comparison with proton/anti- proton 
production data by the TASSO, TOPAZ, DELPHI, and SLD 
collaborations. The scale is Q=32 GeV for theoretical calcu- 
lations in the top figure. 
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B. Optimum fragmentation functions and their 
uncertainties 

The obtained fragmentation functions and their un- 
certainties are shown for 7r+ in Figs. [9] and [TOl The 
gluon, up-quark, and anti-up-quark functions are shown 
at Q^ = l GeV^ in Fig. [9l The charm- and bottom-quark 
functions are shown at their thresholds. The dashed and 
solid curves indicate LO and NLO results, and the dark- 
and light-shaded areas indicate their one-tr uncertainty 
regions estimated by the Hessian method for the LO and 
NLO, respectively. There are differences between the LO 
and NLO functions. The gluon function becomes larger 
in NLO than the LO one, whereas the quark functions are 
smaller in NLO. In the NLO, the favored function D'^^ 
is the largest, and the disfavored one D^^ is smaller than 
D^^ . The gluon function is in-between, and its moment 
{Mg) is roughly given by their average, (M„ + Mu)/2, 
which agrees with the assumption in Ref . [H| . 

Since the experimental data are shown by the sum of 
light-quark flavors, the flavor separation as defined by 
the favored and disfavored initial functions introduces 
uncertainties. For example. Fig. 2] indicates 2% error 
coming directly from the experimental data on the light- 
quark (w, d, s) fragmentation function at 2 = 0.2 and 
Q = Mz, whereas the flavor-separated u-quark function 
in the NLO has 30% error at z = 0.2 in Fig. [TUl In order 
to flnd such a flavor-separation effect on the uncertainties 
at Q = 1 GeV, the data should be fltted by the function 
D'^^ = Nj^^ {l-z)'^'^s where Qs = u + u + d + d+ s + s. 

An error from the assumed functional form is not in- 
cluded in estimating the uncertainty bands. We fixed one 
of the gluon parameters {(3g = 8) , so that the uncertainty 
could be underestimated in zDg^ at large z. However, 

it would not affect the figure of zDg^ as long as /3g ^ 1 
because the distribution itself is small at large z. 

The errors are large in both LO and NLO, which means 
that the fragmentation functions are not well determined 
particularly at small Q^. However, it is important to find 
that all the functions for the pion are determined much 
better in the NLO analysis than the LO ones because the 
uncertainty bands are smaller in Fig. (S] It is especially 
noteworthy that the gluon function is determined well 
in the NLO. The gluonic contributions affect the cross 
section through the NLO coefficient function and NLO 
splitting functions. Therefore, the shrinkage of the er- 
ror band suggests that such gluonic effects are reflected 
in the current inclusive data for the pion. In particu- 
lar, the TASSO collaboration provided many data in the 
small region {Q^ << Aff), and they are important 
for identifying such NLO effects in comparison with other 
data at — M|. The reason why the uncertainties are 
large at small z, especially in the up and charm functions, 
is since small- 2; data are not included in the analysis. 

The functions for the pion are evolved to — 
and the results are shown in Fig. [TOl The functions be- 
come steep ones peaked at z = 0, and the uncertainties 




0.2 0.4 0.6 0.8 1 
z 



0.2 0.4 0.6 0.8 
Z 



FIG. 9: (Color online) Fragmentation functions and their un- 
certainties are shown for at Q^=l GeV^, m^, and ml. 
The dashed and solid curves indicate LO and NLO results, 
and the LO and NLO uncertainties are shown by the dark- 
and hght-shaded bands, respectively. 
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NLO 
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FIG. 10: (Color online) Fragmentation functions and their 
uncertainties are shown for tt^ at — M%. The notations 
are the same as Fig. [5] 
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FIG. 11: (Color online) Fragmentation functions and their 
uncertainties are shown for K'^ at Q^=l GeV^, m^, and m\. 
The notations are the same as Fig. 
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FIG. 12: (Color online) Fragmentation functions and their 
uncertainties are shown for the proton at Q^=l GeV^, m^, 
and ml. The notations are the same as Fig. [9] 



become relatively small in comparison with the ones in 
Fig. [9] However, the gluon uncertainties are still large 
particularly in the LO. The results in Figs. [9] and [TO] 
indicate that the fragmentation functions are not well 
determined in the small region although their uncer- 
tainties are relatively small at large -^^i) especially 
in the NLO. Therefore, if the fragmentation functions are 
used in the small or small pt region << M|) 
such as the lepton-scattering processes \v\ of COMPASS, 
HERMES, and JLab (Thomas Jefferson National Accel- 
erator Facility) and the hadron collisions of RHIC 2] , it 
is very important to indicate the uncertainty ranges of 
the fragmentation functions. 

The obtained fragmentation functions and their uncer- 
tainties are shown for K'^ in Fig. [TT]at Q^=l GeV^, m^, 
and TOj. There are similar tendencies to the pion results. 
The obtained gluon function in the NLO is larger than 
the LO one as the pion case, whereas the quark func- 
tions are smaller. It is interesting to find that the anti- 
strange function is generally larger than the up function: 
D^^ > D^^ , which could be interpreted in the following 
way. In order to create from a parent s (or u), a, uu 
(ss) pair needs to be created. Since the strange-quark 
mass is larger, the ss creation could be suppressed in 
comparison with the uu creation, which leads to the in- 
equality. However, the large uncertainty bands indicate 
that the separation between u and s functions is difficult. 

There is a conspicuous difference between the gluon 
functions for the pion and kaon. The gluon function 



whereas it is at 2 = 0.2 



(zDg ) is peaked at large z, 
0.3 in the pion. Even if an initial distribution with a 
peak at small z is supplied in the fitj the outcome is 
always peaked at large z. It could be physically under- 
stood in the following simple picture. In order to produce 
from a gluon, the gluon should first split into a ss 
pair. Then, another gluon is emitted from the s or s 
quark, and it subsequently splits into a uu pair. It re- 
quires higher energy for the parent gluon to produce the 
ss pair (g ss) in the kaon creation than the one for a 
uu pair {g uu) or dd pair (g — > dd) in the pion creation 
because of the mass difference. The higher energy means 
that the function is peaked at larger z in the kaon. 

The kaon functions also have large uncertainties in 
both favored and disfavored cases. They have slightly 
larger errors than the pionic ones if the ratios SDi/Di 
are considered. The uncertainty bands become smaller 
in NLO than the LO ones. However, the NLO improve- 
ment is not as clear as the pionic one. A possible reason 
is that many accurate data are not taken at small 
(<< M|), for example, by the TASSO collaboration as 
for the pion. 

The fragmentation functions for the proton are shown 
in Fig. [12] at Q^=l GeV^, m^, and mf. Here, the gluon 
moments are fixed by the favored and disfavored mo- 
ments, so that they are almost the same in LO and NLO. 
As expected, the favored functions and are larger 
than the disfavored functions. The gluon functions have 
peaks in the medium-z region. In general, the proton 
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functions are also not determined well, and the uncertain- 
ties are as large as the kaonic ones. The NLO improve- 
ment is also not obvious in the proton. This fact suggests 
that the current proton and antiproton data should not 
be much sensitive to the NLO corrections. 

Since the gluon moment is given by the average of fa- 
vored and disfavored moments, the error of the gluon 
function could be underestimated. A noticeable dif- 
ference from the pion and kaon figures is that the gluon 
uncertainty bands shrink in the region 0.5 < z < 0.7, 
which is caused by fixing the moment. There are large 
contributions from the diagonal terms in Eq. (|39p. The 
term {dDg/dMg)^ is a smooth function of z with a peak 
in the same position as the function Dg. On the other 
hand, (dDg/dag)"^ and {dDg/d(3g)^ have double peaks, 
which come from destructive interferences between two 
derivative terms. For example, the parameter Ng de- 
pends on the parameter Pg by Eq. ([T^ . The term 
dNg/dPg is positive, whereas there is a negative con- 
tribution from the derivative of (1 — z)^!! . In the pion 
case, {dDg/dMg)^H^j^^j^j^ is several times larger than 
{dDg/dag)^H~^^^, so that such double-peak structure 
does not appear in the uncertainty band in Fig. [9l 

^2 ir-i does not exist 

is fixed. 
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However, the term {dDg/dMgYH]^\,j 
in the proton analysis since the parameter Mg 
Therefore, the double-peak shape of {dDg/df3g)'^H 
becomes apparent in Fig. [T^ 

The determined fragmentation functions and their un- 
certainties can be calculated by usin g a library code, 
which is obtained from our web page [3^. The details 
are explained on the web page and distributed files. 

In these analyses for the pion, kaon, and proton, the 
charm and bottom functions are determined mainly from 
the heavy-quark tagging data by the SLD and DELPHI 
collaborations. The light-quark functions are constrained 
especially by the accurate data of the SLD. In order to de- 
termine the fragmentation functions accurately, we need 
low-energy data possibly by the Belle and Babar collab- 
orations. The low-energy data should be important for 
the determination of the gluon functions Dg which are 
essential for describing the low-px hadron productions, 
for example, at RHIC. 

The fragmentation functions for tt^ , K~ , and p can be 
calculated from the obtained functions by using the rela- 
tions in Eqs. ([M]) and ([55)1 . The fragmentation functions 
for 7r°, K^, K^, n, and n can be calculated by using the 
determined functions for the tt"*", K~^, and proton. The 
relations are explained in Appendix \X\ 



C. Comparison with other parametrizations 

Our NLO fragmentation functions are compared with 
other parametrizations by KKP [4], AKK 4], and Kret- 
zer in Figs. [T31 and [15] for the pion, kaon, and 
proton/anti-proton, respectively. Since tt"*", , and pro- 
ton are not separated from tt~ , K~, and anti-proton in 
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FIG. 13: (Color online) The NLO fragmentation functions 
for the pion, z(-D7 + )/2, are compared with other 
parametrizations at Q^=2 GeV^ 10 GeV^, and 100 GeV^ 
The solid, dashed, dash-dotted, and dotted curves indicate 
HKNS, KKP, AKK, and Kretzer parametrizations, respec- 
tively, in the NLO MS scheme. The HKNS uncertainty bands 
are shown by the shaded areas. 
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FIG. 14: (Color online) The NLO fragmentation functions 
for the kaon, z{Df + Df )/2, are compared with other 
parametrizations at Q^=2 GeV^ 10 GeV^, and 100 GeV^. 
The notations are the same as Fig. [TSj 
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FIG. 15: (Color online) The NLO fragmentation functions 
for the proton/anti-proton, + _Df)/2, are compared with 
other parametrizations at Q^=2 GeV^, 10 GeV^, and 100 
GeV^. The notations are the same as Fig. 1131 



the KKP and AKK parametrizations, the combinations 
z{Df + Df )/2, z{Df\ + Df- )/2, and z{Dl + Dl)/2 
are shown. Our parametrization is denoted HKNS (Hirai, 
Kumano, Nagai, Sudoh) in these figures. The hght-quark 
functions are separated in the AKK parametrization due 
to additional OPAL data, which are not used in the other 
analyses. The KKP, AKK, and Kretzer distributions are 
calculated by using their library codes. The small-z part 
is not plotted for the AKK because their code does not 
support the region, z < 0.1045 by considering resum- 
mation effects in comparison with the data. The heavy- 
quark thresholds are taken 4to^ and 4to^ for the fragmen- 
tation functions in the KKP and AKK parametrizations, 
whereas they are and in the Kretzer and HKNS 
analyses. The minimum value of in the KKP and 
AKK codes is GeV^. Therefore, comparisons are 

made at (3^=2, 10, and 100 GeV^. 

If the distributions are compared with each other in 
Fig. [T3]for the pion, they agree well in up-, charm-, and 
bottom-quark functions. The up-quark function of the 
KKP has a singular behavior as z — > 0, which is different 
from the Kretzer's, AKK, and HKNS functions. How- 
ever, the small-z data are not included in these analyses, 
so that the differences should not be taken seriously at 
small z. Since new SLD data in 2004 are accurate and 
they are not used in other parametrizations, the frag- 
mentation functions could be better determined in our 
parametrization. The four types of functions are very 
different in the gluon and disfavored strange-quark func- 
tions in Fig. 1131 All the gluon functions have various 
peak positions in the region 0.1 < z < 0.4 and func- 



tional forms are different. Our gluon function roughly 
agrees with the Kretzer's function. The four strange- 
quark functions are also much different. Although our 
strange-quark function is almost equal to the KKP, the 
Kretzer's (AKK) function is much larger (smaller). How- 
ever, it is important to find that the gluon and strange- 
quark functions, needless to say the up-, charm-, and 
bottom-quark functions, by the KKP, AKK, Kretzer, and 
HKNS are consistent with each other because they are 
within the uncertainty bands. This fact indicates that 
all the analyses are successful and consistent. 

For the kaon, the KKP, AKK, Kretzer's, and HKNS 
results agree well in the charm- and bottom-quark func- 
tions as shown in Fig. 1141 However, they differ much in 
the gluon and light-quark functions. In particular, the 
HKNS and AKK gluon functions are peaked at large z, 
whereas the KKP diverges at small z and the Kretzer's 
hmction is peaked at z 0.4. The KKP and Kretzer 
functions are outside the uncertainty band at small z. 
Since the parameter ag is fixed and because of the func- 
tional form with the peak at large z, the uncertainty 
could be underestimated in the small-z region. The large 
variations among the parametrizations of Dg^ 
indicate that the current kaon data are not accurate 
enough to fix the gluon fragmentation function for the 
kaon from the scaling violation. Lower-energy date, for 
example from Belle and BaBar, should be able to im- 
prove the situation. The light-quark (u, d, s) functions 
are also very different; however, they are roughly within 
the uncertainty band. Except for the gluon function, all 
the analysis results are consistent in the kaon. 

The results agree well with each other in the charm- 
and bottom-quark functions for the proton and anti- 
proton as shown in Fig. 1151 although there are some 
differences at small z where experimental data are not 
used in the analyses. The Kretzer's parametrization is 
not available for the proton. The three gluon functions 
are very different but they are roughly within the un- 
certainty band. It is almost impossible to determine the 
accurate gluon function at this stage. The disfavored 
strange-quark functions are also very different. In partic- 
ular, the KKP (AKK) function is much larger (smaller) 
than our result. Our favored functions, up- and down- 
quark ones, are larger than the KKP and AKK functions. 
However, all these fragmentation functions including the 
disfavored ones are roughly within the estimated uncer- 
tainties. 

From these comparisons, we found that the analyses 
of KKP, AKK, Kretzer, and HKNS are consistent be- 
cause they are generally within the uncertainty bands 
estimated by the Hessian method. However, there are 
noticeable differences in the gluon and light-quark frag- 
mentation functions. They should be clarified by future 
measurements especially at small Q" . The large uncer- 
tainties could cause serious effects in discussing hadron- 
production processes at small pt such as RHIC, HER- 
MES, and JLab. 
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V. SUMMARY 

Unpolarized fragmentation functions for the pion, 
kaon, and proton have been determined in the LO and 
NLO from global analyses of e++e~ h+X data. Their 
uncertainties were estimated by the Hessian method. We 
found the large uncertainties in the fragmentation func- 
tions at small (^1 GeV^) although they become rela- 
tively smaller at high energies {Q^ ~ ^z)- particular, 
the gluon and light-quark fragmentation functions have 
large uncertainties. However, they are determined more 
accurately in the NLO analyses than the LO ones for the 
pion and kaon. Because of the large uncertainties at small 
Q^, it is important that such errors need to be taken into 
account in analyzing the hadron-production data in high- 
energy lepton-proton, proton-proton, and nuclear reac- 
tions. Low-energy e^e~ measurements should be valu- 
able for determining especially the gluon fragmentation 
functions by the scaling violation. A code for calculating 
the determined fragmentation functions can be obtained 
from our web page [s^] . 
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APPENDIX A: FRAGMENTATION FUNCTIONS 

FOR TT°, K°, K°, n, AND n 

The fragmentation functions of 7r+, vr^, K^, K~ , p, 
and p are extracted from the experiment by assuming 
charge symmetry in relating Df^, and Df to the 

corresponding ones Df , Df , and Df by Eqs. ([21 and 
psp . One may need the fragmentation functions for tt*^, 
if", K'^, n, or n in one's studies of high-energy hadron 
reactions. The quark compositions of these hadrons are 



7^°((^i^i-d(^)/^/2), K°{ds), K° {ds) , n (udd) , and n{udd). 
Considering these compositions, we relate the obtained 
functions to the ones for 7r°, K^, K'^ , n, and n [11 . The 
TT*^ functions are given by the averages of 7r+ and 7r~: 

of {z, Q') = l[ of (z, + Df (z, Q2) ]. (Al) 
The functions are also calculated by the relations: 

Df{z,Ql) = Dr{z,Ql)^DT{z,Ql) 
= Df{z,Ql) ^Df{z,Ql), 
Df,iz,Q') = D^:{z^Q^)^D${z,Q% 
Df{z,Q^)^Df{z,Q^). (A2) 

However, one should be careful about the following rela- 
tions at different from Qq because of the NLO evolu- 
tion: 

Df{z,Q') + Df{z,Q') ^ Df{z,Q'), (A3) 

although D^\z,Q'^) = D^°{z,Q'^) is stih valid. Then, 
K'^ functions are related to the K'^ ones by 

Df\z,Q^) = Df{z,Q^), 
Df{z,Q^)^Df{z,Q% (A4) 
The neutron functions are related to the proton ones by 

Dl{z,Q^) = D^,{z,Q^), 
D^iz,Q^) = Dliz,Q^), 

D'l{z,Q^) = DP{z,Q^) fori^u,d. (A5) 

The anti-neutron functions are then given by 

D^{z,Q^)^D^^iz,Q^), 
D^iz,Q')^D^;iz,Q'). (A6) 

Using the code in Ref. [35|, one should be able calculate 
the fragmentation functions for tt"*", tt*^, tt", -ftT"'", 
, , p, p, n, and n at given z and Q^. 
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